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Active transport of sodium and potassium ions by Na,K-ATPase is accompanied by the enzyme 
conformational transition between El and E2 states. ATP and ADP bind to Na,K-ATPase in the El 
conformation with similar affinity but the properties of enzyme in complexes with these nucleotides are 
different. We have studied thermodynamics of Na,K-ATPase binding with adenine nucleotides at different 
temperatures using isothermal titration calorimetry. Our data indicate that p-phosphate is involved in 
complex formation by increasing the affinity of adenine nucleotides to Na,K-ATPase by an order of 
magnitude, while y-phosphate does not affect it. ATP binding to Na,K-ATPase in contrast to ADP binding 
generates a structural transition in the enzyme, which is consistent with the movement of a significant 
portion of the surface area to a solvent-protected state. We propose that ATP binding leads to convergence 
of the nucleotide-binding and phosphorylation domains transferring the enzyme from the "El -open" to 
"El -closed" conformation ready for phosphorylation. 



* These authors 
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I n all animal cells, active transport of sodium and potassium ions across the plasma membrane is facilitated by 
I Na,K-ATPase (sodium pump, EC 3.6.3.9). The transport is followed by conformational transition between the 
I two main states of the enzyme, El and E2, that bind three Na + and two K + ions, respectively 1 . Na,K-ATPase in 
the El conformation has a high affinity to ATP and Na + , while the enzyme in the E2 conformation has a low 
affinity to ATP and high affinity to K + ions 1 . Functional monomer of the Na,K-ATPase consists of the catalytic oc- 
subunit and P-subunit, the latter is required for routing of the a-subunit to the plasma membrane and for the K + 
ions binding 2,3 . The a-subunit contains sites for the binding of nucleotide, Na + and K + , and is responsible for the 
cation translocation as well as for the ATP hydrolysis. A major part of a-subunit is localised in the cytosol and 
consists of the large cytosolic loop forming the nucleotide binding (N) and phosphorylation (P) domains, and the 
small cytosolic loop, which together with the N-terminus of the molecule forms the actuator (A) domain. 
Transition from the El to E2 conformation is accompanied by formation of the phosphorylated intermediate: 
at first y-phosphate of ATP bonded to the N domain is transferred to the carboxyl group of Asp369 located in the 
P domain, and then the acyl-phosphate bond is hydrolyzed 1 . Recently a crystal structure of the protein in 
intermediate [Na 3 ]ElP-ADP state (pig renal a^y enzyme) at 4.3 A 4 and 2.8 A 5 resolution has been solved. 
The structure of Na,K-ATPase in the El free state is absent and it remains unclear which conformational 
transitions the enzyme undergoes when passing from the free El state to the nucleotide-binding state. 

ATP and ADP bind to Na,K-ATPase in the El conformation with similar affinity 6 ' 7 , however ATP and ADP 
affect Na,K-ATPase in a different way. ATP was shown to protect the enzyme from inactivation by hydroxyl 
radicals and ADP does not produce such effect 8 . Structural organization of the free Na,K-ATPase and as com- 
plexes with ATP or ADP may be different. This assumption is supported by the fact that crystal structures of the 
Na,K-ATPase homolog - sarco (endo) plasmic reticulum Ca 2+ -ATPase (SERCA) in the Ca 2+ -bound El state 9,10 , 
and in complex with a nonhydrolyzable ATP analog - adenosine (P - y methylene)-triphosphate n , differ from 
each other by the arrangement of N and A domains with respect to the P domain. 

In this work, we have used isothermal titration calorimetry (ITC) to estimate the equilibrium constants along 
with the enthalpic and entropic components for the A-nucleotides (AMP, ADP and ATP) binding to Na,K- 
ATPase in the El state at different temperatures. Previously the Na,K-ATPase interactions with adenine nucleo- 
tides have been studied by ITC, at a single temperature value and the thermodynamic parameters of the enzyme 
binding with ATP were not reported 6,12 . Thus no conclusions could be made regarding conformational changes in 
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Figure 1 | Na,K- ATPase interaction with adenine nucleotides measured by ITC. Titration curves (upper panels) and binding isotherms (lower panels) 
for duck Na,K-ATPase interaction with AMP (A), ADP (B) and ATP (C) at 37°C and pH 7.5. 



the protein upon binding of nucleotides. In this study we have found 
that P -phosphate is involved in complex formation by increasing the 
affinity of nucleotides to Na,K-ATPase by an order of magnitude, 
while y-phosphate does not affect it. Furthermore, we have demon- 
strated that ATP binding to Na,K-ATPase in contrast to the ADP 
binding generates a structural transition in the enzyme, which is 
compatible with the movement of molecular surface area (ranging 
from 4021 to 5891 A 2 ) from a solvent- exposed to solvent-protected 
state. Comparison of the ITC data with the structure modeling 
results for Na,K-ATPase in the El free state and in the ATP-bound 
state allowed us to conclude that this structural change is caused by 
the locking of the N domain of Na,K-ATPase to the P domain 
induced by the y-phosphate of ATP, thereby reducing the exposed 
surface area of the protein. 

Results 

Affinities of Na,K-ATPase to adenine nucleotides. We used ITC to 
estimate the affinities (association (K a ) and dissociation (K d ) 
equilibrium constants) of AMP, ADP and ATP to duck ATPase at 
different temperatures (Table 1). It should be noted that the 



experimental conditions exclude ATP hydrolysis and formation of 
phosphorylated intermediate (see Methods). A typical set of ITC data 
for AMP, ADP and ATP binding to ATPase at 37°C is shown in 
Figure 1, including the "raw" calorimetric data for the ligand-into- 
protein titration (upper parts) and the binding curves (lower parts). 
The latter were fitted to a model with one adenine nucleotide binding 
site per ATPase molecule. 

In the presence of 3 mM NaCl we observed a 13 -fold and 15 -fold 
higher affinity for the binding of ATPase to ADP compared to AMP 
at 25°C (Ka values equal to 0.07 and 0.9 jiM) and at 37°C (K d values 
equal to 0.04 and 0.6 |iM), respectively (Table 1). Affinities of 
ATPase to ADP and ATP were similar (Table 1). The data show that 
P-phosphate significantly contributes to the affinity of nucleotides to 
ATPase, while y-phosphate of ATP does not affect the affinity during 
complex formation. These data confirm previous suggestions 613 
regarding the key role of P-phosphate in the formation of Na,K- 
ATPase complexes with nucleotides. The obtained binding constants 
for ADP and ATP, are close to the values determined previously 7,14 . 

In the presence of 3 mM KC1 instead of NaCl (E2 state), we failed 
to observe any measurable ATP or ADP binding to ATPase either at 



Table 1 | Thermodynamic parameters of the duck Na,K-ATPase binding to AMP, ADP and ATP determined by isothermal titration calor- 
imetry a 



Ligand 


T, °C 


C M" 1 


Kt. ^M 


AH d , kcal/mole 


TAS e , kcal/mole 


AC?, kcal/mole 


AMP 


25 


1.1 x 10 6 


0.91 


-1.40 


6.83 


-8.23 


AMP 


30 


1.8 x 10 6 


0.56 


-1.81 


6.86 


-8.67 


AMP 


37 


1 .7 x 1 0 6 


0.59 


-2.19 


6.65 


-8.84 


ADP 


25 


1 .4 x 1 0 7 


0.07 


-5.68 


4.06 


-9.74 


ADP 


30 


2.4 x 10 7 


0.04 


-6.00 


4.23 


-10.23 


ADP 


37 


2.6 x 10 7 


0.04 


-6.21 


4.31 


-10.52 


ATP 


10 


1 .2 x 1 0 7 


0.08 


7.70 


16.87 


-9.17 


ATP 


25 


1 .6 x 1 0 7 


0.06 


-4.00 


5.82 


-9.82 


ATP 


37 


2.0 x 10 7 


0.05 


-12.80 


-2.44 


-10.36 



a AII measurements were performed three to four times in imidazole buffer (25 mM ii 

b K a - affinity constant; standard deviation did not exceed ±20%. 

C /Cd - dissociation constant; calculated as 1//C a . 

d AH - enthalpy variation; standard deviation did not exceed ± 1 0%. 

e TAS- entropy variation; calculated from the equation AG = AH - TAS. 

f AG-Gibbs energy; calculated from the equation AG = -RT\r\K a . 



lidazole, 1 mM DTT ; 1 mM EDTA, 3 mM NaCl and 250 mM sucrose, pH 7.5). 
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Figure 2 | Thermodynamic analysis of Na,K- ATPase interactions with 
adenine nucleotides. Temperature dependence of enthalpy of the duck 
Na,K-ATPase interactions with AMP (filled squares), ADP (filled circles) 
and ATP (empty circles) at pH 7.5. 

25 or 37°C, in agreement with the previous data 6 . To determine 
binding constants by ITC it is necessary to use protein concentration 
exceeding the K d value 15 . According to 14 , the K d value for ATP bind- 
ing with Na,K- ATPase in the E2 state is equal to 450 |iM. The con- 
centrations of Na,K-ATPase used in this study render it not possible 
to determine the K d values by ITC. 

Energetics of complex formation between Na,K-ATPase and 
adenine nucleotides. In addition to the binding constant estimates, 
obtained from the shape of the binding curves, we also estimated the 
change of enthalpy, AH, upon binding of an adenine nucleotide to 
Na,K-ATPase (Table 1). Plotting the enthalpy of the interaction 
versus temperature we obtained heat capacity change (AC p ) values 
for the ATPase interactions with adenine nucleotides as d(AH)/dT 
(Figure 2, Table 2). Na,K- ATPase binding both to AMP and ADP 
had low AC p of —66 and —43 cal'mol^'K -1 , respectively, and the 
ATP binding to Na,K- ATPase was characterized by massive AC p of 
— 760 cal-mor^K" 1 (Table 2). On the assumption of additive 
contributions of the different phosphor groups to the AC p values 16 , 
we found that the contribution of y-phosphate is equal to 
-717 cal-mol-^K" 1 . 

Changes in heat capacity are believed to reflect the change in the 
solvent accessible surface area (SAA) during the process of complex 
formation 17 . The large decrease in AC p upon complex formation 
between Na,K-ATPase and ATP suggests that ATP binding to 
Na,K-ATPase significantly decreases its SAA 18 . To assess changes 
in SAA by ITC a number of empirical formulas are used, which 
connect the changes in AC p and the changes of water-accessible polar 
and nonpolar surface area in proteins 19 ' 20 . In the review 21 it was 



Table 2 | Heat capacity changes and change in the solvent-access- 
ible surface area for the duck Na,K-ATPase binding to AMP, ADP 
and ATP a 

Ligand AC p/ cal-mole" 1 -K" 1 AA tot min b , A 2 AA tot max b , A 2 

AMP -66 349 551 

ADP -43 228 333 

ATP -760 4021 5891 

a Heat capacity changes were obtained as d(A/-/)/d7"from Figure 2. 

b Changes in solvent-accessible surface area were estimated using the following formula: AC p = k 
x AA tot/ where AA tot is the total change in solvent accessible surface area in A 2 . AA to t_max and 
AA tot min are calculated for the min and max k values (0. 1 29 and 0. 1 89 cal/(K-mol*A 2 ), as shown 
in 2 ')." 




Figure 3 | Conformational changes of Na,K- ATPase upon ATP binding. 

(A) Cartoon representation of the Na,K-ATPase od-subunit model, built 
on the basis of 2.8 A structure of the porcine od-subunit (Protein Data 
Bank code 3wgu). (B) Complex of the Na,K-ATPase with ATP. The ATP 
was docked to the protein as described 5 . Color coding of the Na,K- ATPase 
refers to individual domains, with the A domain (residues 1 to 77 and 149 
to 270) in yellow, the P domain (residues 363 to 376 and 589 to 753) in 
blue, the N domain (residues 377 to 588) in red, and the transmembrane 
region (M1-M10) in grey. The ATP nucleotide is shown as ball- and- stick 
representation with standard color coding for atom type (carbon, gray; 
oxygen, red; nitrogen, blue; sulfur, yellow; phosphorus, pink; magnesium, 
brown). This figure was prepared with MOE version 2013.08 modeling 
software. 

shown that formulas representing individual contributions of polar 
and nonpolar groups can be replaced by the formula with a single 
coefficient: AC p = k X AA tot , where AA tot is the total change in SAA 
in A 2 . The value of k reported in different studies varies from 0.129 to 
0.189 cal/CK-mol-A 2 ) 21 . Using these values we estimated the total 
decrease in protein SAA when enzyme complexes with nucleotides 
are formed (Table 2). 

Molecular modeling of the Na,K- ATPase structure changes upon 
ATP binding. Models have been constructed of the Na,K- ATPase in 
the "open" El free state and "closed" ATP-bound state (Figure 3). In 
the El "open" state the N and P domains are separated in space. In 
the El "closed" state the N domain is positioned closer to the 
boundary between the P and A domains, with the ATP y- 
phosphate in the N domain interacting with the side chain of 
residue Asp369 in the P domain. The SAAs, estimated for an 
"open" El and "closed" El conformations were 39436 A 2 and 
34687 A 2 , respectively. Thus, the difference in SAA between the 
"open" and "closed" ATP -binding state conformations is 4749 A 2 . 
This is in good agreement with the data of AA tot obtained from ITC 
(Table 2). 

Discussion 

The classical scheme of Na,K- ATPase function consists of the two 
major conformational transitions: between the phosphorylated 
forms of the enzyme E1P-E2P and between the non-phosphorylated 
forms E2-E1. From the crystal structure of Na,K- ATPase it is known 
that in the E2 conformation the N and P domains are positioned 
close to each other 22 , while limited proteolysis results showed that in 
the El conformation these domains reside apart in three-dimen- 
sional space 23 . Conformational transitions in the Na,K-ATPase 
molecule when binding nucleotides up to now has not been 
established. 
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AMP and ADP binding to Na,K-ATPase is followed by small 
changes in the AH and AS values with temperature and there are 
no indications of a structural transition in the enzyme. The large AS 
variation with temperature upon ATP binding to Na,K-ATPase, the 
decrease of AH with increasing of temperature, and the numerically 
large, negative heat capacity change upon complex formation 
(Figure 2, Tables 1,2), suggest that ATP binding to Na,K-ATPase 
induces a structural transition of the enzyme that transfers a part 
of the surface area from a solvent exposed to solvent protected 
state 18 ' 24 . Thus the ATP y-phosphate plays a key role in the structural 
rearrangement essential for the enzyme function. Similarly, binding 
of GTP, but not GDP, leads to a global conformational transition in 
the structures of elongation (EF-G) and translation termination 
(eRF3) factors 19 ' 25 . Modeling of the Na,K-ATPase complex with 
ATP, produces results demonstrating convergence of the ATP y- 
phosphate in the N domain and Asp369 in the P domain of Na, 
K-ATPase (Figure 3B), and shows the change of solvent accessible 
surface area of the protein similar to those that were obtained experi- 
mentally for ATP binding (Table 2). Thus, the structural changes 
occurring upon ATP binding can be attributed to the convergence of 
N and P domains of Na,K-ATPase. For structural homolog of the 
catalytic subunit of Na,K-ATPase, SERCA, it was demonstrated 
that domain arrangements in the Ca 2+ -bound El form and in the 
complex with a nonhydrolyzable ATP analog AMPPCP (Ca2-El- 
AMPPCP) are different 11 . Nucleotide binding assembles the cyto- 
plasmic domains by a 90° rotation of the N domain toward the P 
domain, converting the open Ca 2+ -bound El form into a compact 
structure. 

To summarize, it has been shown that Na,K-ATPase in the El 
form complexed with ATP and ADP have different conformations. 
When complexed with ADP the enzyme does not undergo structural 
changes and remains in the El -open conformation, whereas when 
binding ATP it undergoes a structural transition resulting in a sig- 
nificant decrease of the solvent accessible surface area. Replacement 
of ADP by ATP leads to convergence of the Na,K- ATPase nucleotide 
binding and phosphorylation domains, resulting in the protein 
adopting the El -closed conformation, ready for phosphorylation. 
Thus, y-phosphate is a trigger for the structural change in Na,K- 
ATPase, converting it from the relaxed state to the ratcheted state. 

Methods 

Na,K- ATPase purification and activity measurements. Na,K- ATPase (al(3l 
isozyme 26 ) was purified from duck salt glands (for details see 27 " 29 ) to the purity grade 
of 99% of total protein as confirmed by polyacrylamide gel electrophoresis in the 
presence of SDS. The enzyme preparations also contain a 13 kDa protein that is 
similar to proteins of the FXYD family 30 . Specific activity of the duck Na,K- ATPase 
reached -1800-2000 umol Pi (mg protein X h)" 1 at 37°C. 

Hydrolytic activity of the duck Na,K- ATPase was estimated as ATP cleavage in the 
medium containing (in mM): 130 NaCl, 20 KC1, 3 MgCl 2 , 3 ATP, and 30 imidazole, 
pH 7.4 at 37°C 31 " 33 . Na-K-ATPase-induced ATP hydrolysis was measured as a dif- 
ference in ATP cleavage rate in the presence or absence of the Na-K- ATPase specific 
inhibitor ouabain (1 mM) 29 . Reaction of the ATP cleavage was carried out in the 
volume 0.5 ml. It was initiated by the addition of the protein (2 jag) and stopped 
2 min later by the addition of 0.5 ml of the 3 M cold acetate buffer. The method 
developed by Rathbun and Betlach 33 was used to quantify the amount of inorganic 
phosphate (Pi) spectrophotometrically. In the presence of ouabain, ATPase activity in 
the enzyme preparation was not detected, indicating the absence of other ATPases in 
the preparation. 

Isothermal titration calorimetry (ITC). The thermodynamic parameters of adenine 
nucleotides (ATP, ADP and AMP) binding to duck Na,K-ATPase were measured 
using a MicroCal iTC200 instrument (MicroCal, Northampton, MA), as described 
elsewhere 34 . Experiments were carried out at 10, 25, 30 and 37°C in imidazole buffer 
containing 25 mM imidazole, 1 mM EDTA, 250 mM sucrose, 0.1 mM DTT, and 
either 3 mM NaCl or 3 mM KC1, pH 7.5. Experiments were conducted in the absence 
of Mg 2+ , preventing the formation of phosphorylated intermediate 13 . Absence of ATP 
hydrolysis by the enzyme was monitored by detecting Pi. Aliquots (2.6 ul) of ligands 
were injected into a 0.2-ml cell containing the protein solution to achieve a complete 
binding isotherm. Protein concentration in the cell ranged from 5 to 20 uM, and 
ligand concentration in the syringe ranged from 50 to 200 uM. The heat of dilution 
was measured by injection of the ligand into the buffer solution or by additional 
injections of ligand after saturation; the values obtained were subtracted from the heat 



of reaction to obtain the effective heat of binding. The resulting titration curves were 
fitted using the MicroCal Origin software, assuming one set of binding sites. Affinity 
constants (K a ) and enthalpy variations (AH) were determined and the Gibbs energy 
(AG) and entropy variations (AS) were calculated from the equation: AG = — ^71nK a 
= AH - TAS. 

Modelling. Three-dimensional models of the Na,K- ATPase catalytic a-subunit in the 
El "opened" and "closed" state were created on the basis of the previously published 
2.8 A structure of the porcine a-subunit in the [Na 3 ]ElP-ADP state (PDB code 
3wgu) 5 . For modeling of the complex of Na,K- ATPase in the ATP El "closed" state 
the ATP molecule was placed into the cavity of the ATP -binding site of the enzyme 
similarly to the ADP molecule 5 . Modeling the enzyme in the "opened" El free state 
was carried out by changing the geometry of the polypeptide chain of the boundary 
segments 375-379 and 586-590 of the N-domain, by analogy to the changes observed 
in SERCA upon nucleotide binding, specifically exemplified by the rigid body motion 
of N-domain relative to the A and P-domains 11 . The resulting model structures of Na, 
K- ATPase were minimized in the force field MMFF94x using the MOE version 
2013.08 modeling software (Chemical Computing Group Inc., Montreal, Quebec, 
Canada). The solvent accessible surface area of the a-subunit of Na,K- ATPase was 
calculated with the server http://structure.pitt.edu/anchor. 
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